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Abstract 

O 

' Annihilating dark matter (DM) models based on a scalar hidden sector with Higgs portal-like couplings to 



the Standard Model are considered as a possible explanation for recently observed cosmic ray excesses. Two 
versions of the model are studied, one with non-thermal DM as the origin of the boost factor and one with 
Sommerfeld enhancement. In the case of non-thermal DM, four hidden sector scalars which transform under 



■ a U(l)x symmetry are added. The heaviest scalars decouple and later decay to DM scalars, so providing 

the boost factor necessary to explain the present DM annihilation rate. The mass of the annihilating scalars 

D ; 

is limited to < 600 GeV for the model to remain perturbative. U(l)x breaking to Z2 at the electroweak 



transition mixes light 0(100) MeV hidden sector scalars with the Higgs. The DM scalars annihilate to these 
light scalars, which subsequently decay to two pairs via Higgs mixing, so generating a positron excess 
without antiprotons. Decay to n + n rather than e + e~ is necessary to ensure a fast enough light scalar decay 
rate to evade light scalar domination at nucleosynthesis. In the version with Sommerfeld enhancement 
only three new scalars are necessary. TeV scale DM masses can be accomodated, allowing both the higher 
energy electron plus positron excess and the lower energy PAMELA positron excess to be explained. DM 



annihilates to 2/j 1 /u pairs as in the non-thermal model. This annihilation mode may be favoured by recent 
observations of the electron plus positron excess by FERMI and HESS. 
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I. INTRODUCTION 



Recent results from the satellite experiment PAMELA yj] indicate an excess of positrons at 
10-100 GeV as compared with the expected galactic background, confirming the earlier results 
from HEAT Q2Q and AMS yfl. Surprisingly, PAMELA did not find any antiproton excess below 
100 GeV as compared with the galactic background [4]. Evidence was obtained from the balloon 
experiments ATIC fl] and PPB-BETS fl] of an excess electron plus positron flux as compared 
with the galactic background in the energy range 100-800 GeV. These results have recently been 
reconsidered by FERMI [7] and HESS [8], which do not confirm the large excess and spectral 
features observed by ATIC and PPB-BETS. However, HESS does not rule out the possibility of 
an electron plus positron excess, although there is no indication of structure in the electron plus 
positron spectrum [8fl, while FERMI observes a flattening of the electron spectrum relative to 
that predicted by a conventional diffusive model for the background, which suggests new physics, 
although again no prominent spectral features are observed [7]. Therefore an electron plus positron 
excess remains a possibility. These results raise the exciting prospect that the positron and the 
electron plus positron excesses could be attributed to annihilation of dark matter (DM) particles 1 . 
If DM annihilation is the explanation for the positron excess at lower energies and the possible 
electron plus positron excess at higher energies, then the annihilation rate of DM at present should 
be larger than that expected from the canonical thermal relic annihilation cross-section in the case 
of a smooth distribution of DM in the galaxy (^3x 10 cm 3 /s) II 1411 . This is the boost factor 2 . 
The origin of the boost factor could be astrophysical, because of the merger of sub-structures, or 
entirely from particle physics, or a combination of the two. 

A popular method to achieve the boost factor is Sommerfeld enhancement of the DM anni- 
hilation cross-section Ill4l - ll7ll . This typically requires the introduction of new light bosons of 
mass Mb ~ CuMom m order to mediate a force between the DM particles, where Mom is the DM 
particle mass and a is the interaction's fine-structure constant. (An exception is discussed in [15], 
where the enhancement is mediated by electroweak interactions.) An alternative approach is to 
use non-thermal production of DM to accomodate a large annihilation cross-section, usually via 



1 Nearby astrophysical sources 1191 41 1H and decaying DM II121I13H are also possible explanations. 

2 Different authors define the boost factor in different ways, with some reserving this term for the astrophysical boost 
due to clumpy DM. We will use it to refer to the total enhancement of the DM annihilation rate. 
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heavy particle decay 3 fl 1 8M . In the following we will consider both possibilities 4 . 

In addition, in order to produce positrons without a sizable amount of antiprotons, a mechanism 
to allow DM to annihilate primarily to leptons is required. One approach is to introduce new 
'leptophilic' couplings of the DM particles to leptons uM- In this paper we will instead consider 
a DM sector which interacts with the SM via generic (non-leptophilic) couplings. Our goal is to 
determine whether such non-leptophilic models can account for the observed cosmic ray excesses 
and to obtain the necessary conditions on their masses, couplings and field content. Our analysis 
of the ingredients required to construct successful non-leptophilic models may then guide the 
construction of more complete models which can explain the necessary features. 

Our models are based on a SM singlet scalar sector interacting with the SM via Higgs portal- 
like interactions. Adding a scalar DM particle 5 is a particularly simple way to extend the SM to 



account for DM [24, 



2511 . The scalar is typically stabilised by either a discrete Z2 or U(l) symme- 



try. It interacts with the SM sector via the coupling S^SH^H (which has come to be known as the 



'Higgs portal' D26D), which is the only renormalizable coupling of the S to SM particles. Several 
DM models based on this type of coupling have been proposed [ 2v|] . However, this coupling alone 
cannot account for DM annihilation primarily to leptons, nor can it account for the boost factor. 
Here we extend the symmetry of the DM particle to a sector of SM singlet scalar fields. The DM 
sector is composed of the DM scalar S and additional scalar fields % ; , all of which carry non-trivial 
charges under a symmetry U(l)x- In the model with non-thermal DM, three scalars % ; , i = 1, 2, 3, 
are introduced. The heaviest scalar Xi populates the number density of DM, so providing the boost 
factor, while the lightest scalar %3 ensures an annihilation channel of DM to two ^ + /u pairs. The 
role of %2 is to mix %3 with the Higgs via its vacuum expectation value (VEV), inducing its decay 
to leptonic final states. In order to ensure that %3 decays before dominating at nucleosynthesis, 
U(l)x must be broken to a discrete Z2 symmetry which maintains the stability of S. This occurs 
spontaneously at the electroweak phase transition, when %2 acquires a VEV. U(l)x breaking also 
mixes the lightest scalar %3 with the Higgs boson, providing a mechanism for leptonic annihilation 
of 5 to 2 fj + V pairs 5 . In this version of the model the mass of S is constrained by perturbativity to 
be less than approximately 600 GeV In the version of the model with Sommerfeld enhancement, 
the light scalar %3 which mixes with the Higgs is also used to mediate the enhancement. This 



i We note that in SUSY models a natural alternative is Q-ball decay 

4 Another possibility, annihilation close to a pole, has been considered in 12 ll 12211 . 

5 The Higgs mixing mechanism was first described in 11611 . 
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model requires only two scalars, %2 and %3, and can accomodate a larger range of S mass. DM 
again annihilates to 2 pairs via Higgs mixing as in the model with non-thermal DM. An- 



290, with 



nihilation to muons appears to be favoured by recent data from FERMI and HESS [28, 
annihilation to 4/j being favoured by the analysis of |290. 

Our paper is organised as follows. In Section II we present and discuss a Higgs portal model 
with non-thermal production of DM as the source of the boost factor. In Section III we discuss 
a version of the model with Sommerfeld enhancement in place of non-thermal production. In 
Section IV we present our conclusions. 



H. A HIGGS PORTAL MODEL WITH NON-THERMAL DARK MATTER 



A. Overview 



We extend the SM by adding a dark sector composed of a singlet DM scalar S of mass Ms and 
three additional scalar fields i = 1,2,3 with masses m ; , i = 1,2,3, such that m\ 3> 3> m^. 
We impose a symmetry U(l)x, under which the fields % ; carry a charge +1 and S carries a charge 
+3/2. In order to avoid a Goldstone boson from U(l)x breaking we will consider this to be a 
gauge symmetry. The gauge interaction will not have any significant effect on the cosmological 
evolution of the model, only contributing to the already rapid annihilation and scattering between 
the hidden sector scalars. The SM fields are neutral under U(l)x- The dark sector fields interact 
with the SM via Higgs portal-type couplings to the Higgs bilinear H^H. U(l)x will be broken at 
the electroweak (EW) phase transition to a surviving Z2 symmetry under which S is odd while rest 
of the fields, including the SM fields, are even. Since S is odd under the surviving Z2 symmetry, it 
is stable and a candidate for DM. 

Our model is based on generic couplings of the gauge singlet scalars to H*H. The renormaliz- 
able couplings of the scalar sector of the Lagrangian are given by 

L 2 m 2 ahi+M 2 s S^S + m 2 H H^H 

+ X s (S*S) 2 + X H (H Ji H) 2 +yS f SH*H 

+ T\ijuxhjxhi + (ns)ifXdXjS f s + (n\ H )ijxhj H ^ H > 0) 

where H is the SM Higgs and i = 1,2,3. We assume that all couplings are real and that all particles 
in the dark sector have positive masses squared (m?,M| > 0). As usual mjj < so that H acquires 
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a vacuum expectation value (VEV). If (t\h)22 < and is sufficiently small then %2 gains a 
U(l)x -breaking VEV when the electroweak (EW) phase transition occurs. 

The cosmological evolution of the model can be summarized as follows. %\ is assumed to have 
the largest mass in the hidden sector. S has a mass > 100 GeV if its annihilation is to account 
for the PAMELA observations and so freezes-out at a temperature > 5 GeV. A key requirement 
of the model is that X\ decay occurs sufficiently long after the S density has frozen out of thermal 
equilibrium that it can boost the S relic density. %\ can decay to x]XjXh Xi^S or XiH^H (i, j, k ^ 
1). After the EW phase transition, when Xi acquires a VEV, Xi a l so decays through the two body 
processes: Xi XjXk, with j,k^\ and X\ — >■ S^S, hh, Xjh with j ^ 1 , where h is the physical Higgs 
scalar. We will see that late decay of Xi requires these couplings to be very small, of order 10~ 10 . 
Due to t/(l)x-breaking, X3 mixes with the physical Higgs scalar h and decays to SM fermions via 
the Yukawa couplings. If the X3 mass is in the range 2ra„ to 2m n o (212-270 MeV) then it decays 
predominantly to pairs. For larger X3 mass, X3 decay to n° pairs produces photons while 

%3 decay to nucleon-antinucleon pairs produces antiprotons. Whether the photon flux from pion 
decay is excluded depends on the nature of the DM halo, with cuspy NFW halos excluded but 
cored isothermal halos still likely to be consistent with present bounds 6 [3(1 31]. In the following 



we will consider the X3 mass to be in the range Im^ to 2m 7I o, although the upper bound may be 
increased to 2m pmton if the photon flux from pion decay is within observational limits. This range 
of %3 mass requires that the couplings of X3 to H and %2 are less than 10 . The lifetime for %3 
decay to is short compared with the time at which nucleosynthesis begins, so the relic X3 

density, which would otherwise dominate at nucleosynthesis, safely decays away. This would not 
be true for decay to e + e~ , which would apply if the %3 mass was less than 2m jU . The present S 
density annihilates primarily to X3 pairs which promptly decay to muons. The subsequent decay 
fj + — > e + + v e + v M then accounts for the positron flux without any antiproton flux. 



6 A possible problem with photons from pion decay was previously noted in the context of an axion decay model for 
cosmic ray anomalies 11311 . 
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B. Electroweak Phase Transition and Spontaneous Breaking of t/(l)x 



After the EW phase transition H develops a VEV, which triggers a VEV for % 2 . The VEV of H 
and %2 also induce a VEV for %i through the couplings TI1222 and (t]h)i2 

(xi)Rj _ T ll222 ^ + (T1g ) 12 „ v2 ^ 

where (% 2 ) = w and (//) = v. As we will show, T) 1222 and (r)#) 12 are required to be no larger than 
0(1(T 10 ) in order to ensure the late decay of %\. Therefore with u ~ v ~ 100 GeV and m\ ~ 0(1) 
TeV we find that is negligibly small, 0(100) eV Similarly, %3 also gains a VEV 

, v _ T12333 U 3 + ("l) H )23 UV 2 
m 3 

We will show later that the mass of the lightest mass eigenstate %' 3 must be less than 0(1) GeV 
in order to ensure that it will decay primarily to leptons. This is most easily understood if all the 
terms in the mass matrix are less than 1 GeV, which in turn requires that all the couplings of %3 
to %2 and H are less than O(10~ 6 ). (Larger entries in the mass matrix are possible but would 
require sufficient cancellation between the contributions to the lightest mass eigenvalue.) Since 
m3 is also no larger than 0(1) GeV, this means that (5^3) < w, v. Although a value for < %3 > 
which is comparable to (H) and (% 2 ) is possible, this will not qualitatively alter our results from 
the case where (^3) w, v, since it will only alter the admixtures of % 2 and H in the lightest 
mass eigenstate by 0(1) factors. Therefore, to simplify the analysis we will set < %3 >= in the 
following and consider 

<Xi> = 0, (5)=0, (H)=v, (X2>=« (X3>=0, (4) 

with u and v obtained by minimizing the scalar potential 

V = m 2 H v 2 + mju 2 + r|2222« 4 + k H v 4 + {y\ H )i2U 2 v 2 . (5) 

We assume that (t|//) 22 and are negative with all other terms positive. Vacuum stability re- 
quires that (r|#)22 > — 2a/ti2222^//- Minimizing Eq. © gives 



(2l H m 2 -^ H ) 22 m 2 ) 

(( r l//)| 2 — 4r l2222^//) 
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and 

' (2ri2222W^ - (y\ H ) 2 2m\) 



(7) 



((ri//)^ 2 -4ri2222^H) 

In the following we will assume that r|2222, |Cn//)22| are ~ 0.1 and that m,2 H 100 GeV, therefore 
u ~ v. 

The %2 expectation value breaks the U(\)x symmetry to a discrete symmetry, under which the 
scalars transform as <|> — > e l2% ®§, where Q is the U(l)x charge of (|). Thus %i (2=1) transforms 
as %i —> Xu while S (Q = 3/2) transforms as S — > —S. Therefore 5 is stable due to the residual 
discrete symmetry. 

Once %2 gains a t/(l)x-breaking expectation value, the Higgs h (defined by Re(H°) = (h + 
v) /y/2) mixes with the real parts of %\, %2 and % 3 . The % 2 expectation value is assumed to be ~ v, 
so %2 — h mixing will be large, which should have significant consequences for Higgs phenomenol- 
ogy. %1,-h mixing provides the mechanism for DM to annihilate primarily to lepton final states. 
For simplicity we will consider only the mixing of the Higgs with %3, which is responsible for the 
important physics. In the basis spanned by VlRe%3 and h, the effective mass squared matrix is 
given by 

/m^ + 6T|2233" 2 + (Tl//)33V 2 (T[h)2?>UV \ 

M 2 =\ . (8) 

V 0lff)23wv 2X H v 2 + (r[ H ) 2 2U 2 J 

In this we have assumed that r^z is independent of the order of i, j, k, I. Diagonalising this gives 
mass eigenstates % 3 and h with 

Mj « m 2 + eTto, 2 + (^) 33 v 2 - iM^l (9 ) 

h 

and 

M 2 hf ^2l H v 2 + (T} H ) 22 u 2 , (10) 
where we assume %3 — /z mixing is small. The %3 — /z mixing angle is 

P-^2— • (ID 

h 

In order to have DM annihilation to muon final states we require that M i is in the range 212-270 

X3 

MeV. Therefore we require that m 3 < 300 MeV, T|2233 < 10~ 6 , (11/7)33 < 10~ 6 and {x\ H )23 < 10~ 6 , 
assuming that u ~ v ~ ~ 100 GeV. Although for simplicity we considered only the mixing 
between % 3 and h, in general % 3 will mix with %2 and %\ in addition to % 3 . The corresponding 
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couplings will also be constrained by the requirement that the light eigenstate mass is ~ 0(100) 
MeV, therefore the additional mixings will not change the model qualitatively. This illustrates an 
important feature of generic Higgs portal models for cosmic ray excesses: some couplings must 
be strongly suppressed. We will show that suppressed couplings are also necessary to produce the 
boost factor via non-thermal production of DM. 



C. Non-Thermal Production of 5 Dark Matter 



The S density is due to out-of-equilibrium decay of %\. This must occur at a sufficiently low 
temperature that the S scalars can have a boosted annihilation cross-section without annihilating 
away after being produced by %i decay. An initial density produced by %\ decay at temperature 
Tdecay will annihilate down to a density at T^ ecay given by 



„ H(Tdecax i . , 

ns{Tdecay) = 7-j fT" , (12) 

(0|Vrel|)5 



where ns is the number density and (o|v re i|)s is the annihilation cross-section times relative ve- 
locity (which is T independent for the case of annihilating scalars). Eq. (fT2l) is true if the initial S 
number density from %i decay is larger than nsijdecay)- Since ns g(T)T 3 while H <* g(T) l / 2 T 2 , 
where g(T) is the effective number of relativistic degrees of freedom, Eq. (fT2]) implies that the 
present ratio of the S number density from %i decay to that from thermal freeze-out (which is 
given by Eq. (fl2b with Tdecay replaced by the S freeze-out temperature T$) is 



nsdecay _ ( g(T S ) \ 1/2 T S 



(13) 



n Sth \g (Tdecay) J Tdecay 

Since the annihilation cross-section is enhanced by a factor B over that which accounts for ob- 
served DM via a thermal relic density, it follows that the thermal relic S density is smaller than 
that observed by a factor B. So in order to account for the observed DM via %\ decay we must 
require that n Sde cay ~ B xn Sth . Therefore 

( ggg) y /2 Ts 



TdeCay ~{g(Tdecay)) B~ ^ 

In this we have neglected the logarithmic dependence of the freeze-out temperature on the annihi- 
lation cross-section and therefore treated T$ as a constant. Ts is related to the S mass by T$ = M$/zs 
with zs ~ 20 [33]. For example, if ms ~ 400 GeV then Ts ~ 20 GeV. Since the observed positron 
and electron excess requires that B ~ 10 3 , we would then require that Tdecay ~ 20 MeV. 
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This very low decay temperature is difficult to achieve via particle decay, as it implies a very 
long-lived particle with lifetime x ~ H 1 ~ lO -3 ^. %i can decay to S pairs via the three-body 
decays %i — > %2S^S and % 3 S^S. The decay rate is given by 

r *. 53 m? M t • (15) 

where X 2 = (v\s) 2 2 + ( r l5)i3- I n order to have a late enough %i decay we require that T Xl < H at 
T = O(10) MeV, which implies that X < 10~ 10 . %i can also decay to S pairs via the two-body 
decay %\ — > S^S once %2 gains a VEV. The decay rate is given by 

r w Mij^ . (16) 

With w ~ 100 GeV and ~ 1 TeV, this decay rate as a function of (115)12 is comparable to 
Eq.(fT5l) as a function of X. Therefore we also require that (y\s)n ~ 10~ 10 . (Similar constraints 
apply to other %i two-body decay modes.) 

The above assumes that the initial S density from %i decay is larger than that given in Eq. (PT21) . 
This requires that n %l {T decay ) > n s (T decay )/2, since each Xl decay produces 2 S. n %l (T decay ) is 
given by 

frp \ _ g(T deca y) T decay ^(^%i) (\H\ 
n %\ \ L decay) ~ 7^~\ ~3 TZi |\ ) V 1 ' ) 

siW Tl (a|v rel |) Xl 

where T Xl is the Xl freeze-out temperature. The condition n Xl (T aecay ) > ns(T decay )/2 then implies 
that 

g{T decay ) 1 / 2 T decay 1 (g|Vrel|)xi ^n, 

g(T Xi y/ 2 T Xl 2 (o|v re i|)s ' 
This translates into an upper bound on B, 

8 <« 15 «^^. (19) 
(CM)*, zs M x , g(T Xt )m 

Therefore if B < Bq ps (cj|v re i|)s/ (^KeiDxi (assuming that M %1 ~ M5 and z Xl ~ zs) then the re- 
quired %i decay temperature will be given by Eq. (fl4)) . The cross-section times relative velocity 
for non-relativistic Xl pair annihilation to S and // is given by 

( ° |Vrel|)xi = 32^M[ ^ s)2n + {J}H)2 ^ ■ (20) 
The annihilation cross-section times relative velocity for non-relativistic S is given by 

( ° kel|)5= 3i^ + 72] ' / = 2 ' 3 • (21) 

s 



With B ~ 10 2 - 10 3 , Eq. (OH) is therefore satisfied if (1^)11,(^)11 < 10" 2 , assuming that y ~ 0.1. 
If this is not satisfied and B > Bq, then the boost factor is given by Bq rather than B. In this case the 
S density comes directly from %\ decay without subsequent annihilations. An even lower %\ decay 
temperature and smaller (^5)12 would then be necessary in order to account for the observed DM 
density. 



D. Dark Matter Annihilation Rate 

The dominant S annihilation mode is assumed to be to %' 3 pairs. In this case the annihilation 
cross-section times relative velocity is given by 



32tcM 2 



^hVe,|),= T^|- (22) 



In order to account for the cosmic ray excesses, the annihilation cross-section times relative veloc- 
ity necessary to account for thermal relic DM, (o|v re /|) 3 x 10~ 26 cm 3 /s = 2.6 x 10~ 9 GeV 2 , 
must be boosted by B ~ 10 2 - 10 3 for DM masses in the range O(100)GeV - 0(l)TeV. This re- 
quires that 

(tis) 3 3«(5-16)x (V^) • (23) 



\ TeV, 

Therefore for the theory to remain perturbative ((115)33 < 3) we require that Ms ^ 600 (190) GeV 
for B = 10 2 (10 3 ). Thus while the model can account for the positron excess in the range 1-100 
GeV observed by PAMELA, an electron plus positron excess at energies up to 0(l)TeV cannot be 
explained if the model is to remain perturbative, in which case an alternative explanation for the 
electron plus positron excess is required, most likely astrophysical. This conclusion is likely to 
apply rather generally to models which do not have Sommerfeld enhancement of the annihilation 
cross-section. (However, this does not exclude the possibility of a large annihilation cross-section 



due to strong coupling between S and %3, which is only constrained by unitarity [32]].) 



E. Leptonic Final States via Dark Matter Annihilation to % 3 

In order to account for the positron excess without an accompanying antiproton flux, the S 
annihilations at present should proceed primarily through leptonic decay channels. In our model 
this is achieved through a mixture of %3-Higgs mixing and kinematics. U(l)x breaking due to 
< %2 > causes the real part of %3 to mix with h. If the dominant S annihilation mode is S^S — > %3 t %' 3 
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and if the mass is in the range 2m u < M > < 2m % o, then the mixing of %3 with h leads to the 

X3 

decay % 3 — > via the muon Yukawa coupling. This is illustrated in Figure 1, treating the 

mixing as a mass insertion. Thus S annihilation will produce a A[x final state via the process shown 
in Figure 2. 

The decay rate for % 3 — > is given by 

r,=^, s (D5Mf»y«„ (24) 

x 3 871 x 3 871 \Mjy fe' 
where is the Yukawa coupling of the SM Higgs to fJ + fJ ■ This gives for the lifetime of %' 3 

X/ = — = 4 x 10 4 - — — 

*3 T x , \(t\h)23J V ^ 



/ M h \ 4 / 200Mev\ 



where we have used m = 100 GeV and v = 174 GeV. The short lifetime of %' 3 ensures that the ther- 
mally produced %' 3 will decay well before the onset of nucleosynthesis at 0(l)s. This is essential, 
as the %' 3 will freeze-out while relativistic (since there are no annihilation channels for %' 3 once 
T < M X2 ,M,/) and so they will dominate the energy density at nucleosynthesis. The decay rate to 
e + e~ is suppressed by Y^/Y^ ~ 10~ 5 , leading to a lifetime ~ 10 s. Thus for %' 3 to decay before 
nucleosynthesis, decay to muon pairs must be kinematically allowed 1 . 



<*2> 




FIG. 1: %3 decay into /h + ili via %3 -Higgs mixing and the muon Yukawa coupling. 



7 We note that the mechanism described in [16], which is based on a single scalar field (j) and which gives a deca y ra te 
equivalent to that here but with < X2 > replaced by < %3 >^ 200 MeV (where 5(3 is the equivalent of (j) in 11611 ). 
results in a lifetime which is too long and so (j) domination at nucleosynthesis. 
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FIG. 2: The primary S annihilation process to 4/j via %' 3 decay. 

The previous discussion applies to the decay of the real part of %3, which mixes with the physi- 
cal Higgs. The imaginary part of %3 does not mix with the Higgs and is therefore stable. However, 
so long as the %3 self-coupling ^3333 is large, the %3 scalars will maintain a thermal equilibrium 
with each other even after they have decoupled from thermal equilibrium with other particles. 
Since the imaginary part of %3 is generally heavier than the real part after the latter mixes with the 
Higgs, the imaginary part of %3 will annihilate to the lighter real part once T < M X3 , thus ensuring 
that the entire %3 density can decay to muons prior to nucleosynthesis. 

F. Sub-dominant 5 annihilation to Higgs pairs 

We have so far considered the annihilation S^S — >■ via the quartic coupling (1^5)33, which 
primarily produces pairs. However, it is also possible to have S^S — > H^H via the cou- 

pling 8 y. The branching ratio to Higgs pairs B s ^ s ^ H i H w WCns)^ should be small enough 
that the production of Higgs pairs does not result in a large antiproton signal. This requires that 
B S f S ^ H t H ^0.1. However, this still allows a significant coupling to Higgs pairs, which can con- 
tribute a small antiproton component to the cosmic rays from DM annihilation. The S^SH^H 

8 We are assuming that the S mass is sufficiently large that we can approximately calculate the annihilation rate in 
the < H limit. 
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coupling also mediates the coupling of S to nucleons, which may allow direct detection of DM. 
These possibilities are distinctive features of Higgs portal models which distinguish them from 
models based on purely leptophilic couplings. 



G. Positron Excesses from S^S Annihilation 



The annihilation of S^S pairs will give rise to mostly /j ± pairs which finally decay to e ± and 
neutrinos. The electrons and positrons from S^S annihilation then travel under the influence of the 
galactic magnetic field and therefore the motion of e ± is expected to be a random walk. As a result 
a fraction of e ± flux will reach the solar system. 

The positron flux in the vicinity of the solar system can be obtained by solving the diffusion 



equation [[15 



34 



35] 



jjAEJ,t)=KAE)V 2 fe+(E 1 r,t) + j t [b(E)fAEJ,t)]+Q{E 1 r) 1 



(26) 



where f e +(E,r,t) is the number density of positrons per unit energy E, K e +(E) is the diffusion 
constant, b{E) is the energy-loss rate and Q(E,r) is the positron source term. The positron source 
term Q(E,r) from S^S annihilation is given by 

t dN P + 



Q(E 1 r)=n 2 s (r)o s \v K i 



dE 



(27) 



In the above equation the fragmentation function dN e +/dE represents the number of positrons 
with energy E which are produced from the annihilation of S^S. 

We assume that the positrons are in a steady state, i.e. df e + /dt = 0. Then from Eq. (1261) . the 



positron flux in the vicinity of the solar system can be given in a semi-analytical form [[15 
*AEf®) = 4^te)!<*M J E dE'-^flCK D (E,E')) , 



34, 



35] 



(28) 



4%b(E) 

where Xd{E,E') is the diffusion length from energy E' to energy E and I(Kd{E,E')) is the halo 
function which is independent of particle physics. An analogous solution for electron flux can also 
be obtained. 

Positrons in our galaxy are not only produced by S^S annihilation but also by the scattering 
of cosmic -ray protons with the interstellar medium [36]. Thus the positrons produced from the 
latter sources can act as background for the positrons produced from the annihilation of S^S. 
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The background fluxes Bol of primary and secondary electrons and secondary positrons can be 
parameterized as 0711 : 

0.16e 1 1 



i 

3> 



bk g 



prim, e 
bkg 



bkg 



l + lle°- 9 + 3.2£ 2 -i5 

0.70E - 7 
1 + 1 le 1 - 5 + 600e 2 - 9 + 580e 4 2 
4.5£ a7 



GeV^cnrV 1 !^ 1 



GeV~ 1 cm~ 2 s~ 1 sr~ 1 



GeV~ 1 cm~ 2 s~ 1 sr~ 1 



(29) 



l+650£ 2 - 3 + 1500e 4 - 2 

where the dimensionless parameter £=E/(1 GeV). The net positron flux in the galactic medium is 
then given by 

(^ e +) G ai = (<V)bk g + <V(£,F ) • (30) 

The second term in the above equation is given by Eq. (l28l) . which depends on various factors: 
b(E), \d{E,E'), I{X d {E,E')), v e +, (hs)q and the injection spectrum dN e +/dE' . The energy loss 
term b(E) (due to inverse Compton scattering and synchrotron radiation due to the galactic mag- 
netic field) is determined by the photon density and the strength of magnetic field. Its value is 
taken to be b(E) = 10 _16 £ 2 GeVs _1 |37|1 . The number density of S DM in the solar system is 
given by 

(» S )o = f§, (3D 



where p ~ 0.3GeV/cm 3 . In the energy range we are interested in, the value of v e + is taken 
approximately to be c, the velocity of light. The values of diffusion length ^dULiE!) and the 



corresponding halo function I{X D {E,E')) are based on astrophysical assumptions II 15l. l34l 13511 . By 



considering different heights of the galactic plane and different DM halo profiles the results may 



vary slightly. In the following we take the height of the galactic plane to be <, 4 kpc, w 



lich is 



referred to as the "MED" model [[15l 



34 



35|], and we have used the NFW DM halo profile [38], 

2 




(32) 



to determine the halo function I(Xj)(E,E')), where r s ps 20kpc and r ^ 8.5kpc. (We find that our 
results are not strongly sensitive to the halo profile.) In Figure [31 plotted using DARKSUSY [39], 
the positron fraction from S^S annihilation is compared with the data from AMS, HEAT and 
PAMELA for the case of Ms = 600 GeV, showing that a good fit is obtained in this case. 
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FIG. 3: Positron fraction from S^S — > 2^ + n at Ms = 600 GeV. We have used the annihilation cross-section 
(a|v|) =4.5 x 10" 23 cm 3 /s. 



H. Nucleosynthesis, Gamma-Ray and CMB Constraints on Enhanced S^S Annihilation 

So far we have considered large annihilation cross sections of the order of 10~ 23 cm 3 s _1 in 
order to fit the excess of the observed cosmic-ray electron fraction. This value is approximately 
10 — 10 3 times larger than the canonical value of the annihilation cross section for thermal relic 
DM (~ 3 x 10~ 26 cm 3 s _1 ). Therefore we have to check if this value is consistent with other 
cosmological and astrophysical constraints, in particular those from nucleosynthesis and due to 
gamma-rays from the galactic centre (GC) and halo. We will consider S annihilation primarily 
to %' 3 pairs, but we will include the possibility of a small but significant branching ratio to Higgs 
pairs. 



First we shall discuss constraints which come from BBN [40-43]. Even after the freeze-out 
of S^S annihilations, a small amount of S pairs continue to annihilate. In our model, the S^S pair 
dominantly annihilates into pairs with some fraction into H^H. 

The photodissociation of D and 4 He is severely constrained by observational value of 3 He/D. 
According to |42j] we have a constraint on the annihilation cross section into pairs, 



/ \ W \photo 1-0 xlQ- 21 3 _ l f E vis /M s \ V M s \ 
(ov) < (cv>p , , = cm 3 s — — — — , (33) 



'sr^VAr Bfs^r V 0.7 ; VlTeV 

where E v - ls /Ms represents the fraction of the total energy 2Ms which goes into visible energy E V - 1S 
i.e. charged particles and photons. #sts_>.2ju+ir ^ s me branching ratio for S^S annihilation into 2 
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FIG. 4: Upper bounds from BBN on the annihilation cross section of S'S into a Higgs H'H pair as a 
function of the DM mass, where the branching ratio is normalized to B s t s ^ H t H = 1. Here we have assumed 
the mass of Higgs boson is 130 GeV. The name of the light element used for the constraint is written near 
each line. The vertical band at the left side indicates the region which is not kinematically allowed. 



pairs. In case of the muon decay, E y [ s ~ 0.7Ms- 
In addition, in order to limit the branching ratio to Higgs pairs, we have calculated the con- 
straints on the cross section to H^H which follow from photodissociation and hadron emission. 
For photodissociation we find [|42|l 

/ \w \Photo 7.0 xl(T 22 3 _j fE wis /M s \ _1 ( M S , 



where PYTHIA [44] gives 

^Vis ~ 1.0A/<5 and and Bgf$_ ^jjfjj is the branching ratio into H^H. (In 



the low energy limit this becomes the total branching ratio to W, Z and h.) In the current case 
B S i s-^h^h = ^ — ^s 1 s^2/j+fj- ■ The dominant upper bound comes from the smaller of (^^^^fj+fj- 
and (^^t^z/iT/- These bounds are generally compatible with the range of values required to 
account for the cosmic ray excesses, (ov) ~ 3 x 10~ 24 — 3 x 10 cm s for boost factors 10 2 — 
10 3 . 

The most severe bound on hadron emission comes from the overproduction of deuterium by the 
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destruction of 4 He. This process is constrained by observational D/H Il42h . From this we obtain 
the hadron emission constraint 

/ v / vhad 1.3xl(T 23 3 

with N n the number of emitted neutrons per single annihilation. In the case of H^H emission, N n is 



approximately 1.0, which is obtained using PYTHIA Q44Q . This is again consistent with the range 
of (cv) required to account for the observed cosmic ray excesses. We have plotted the results for 
annihilation into H^H in Figure 4, with the normalization B S \ S ^ H ^ H = 1. 

In summary, for the range of S mass which is compatible with perturbative couplings, the boost 
factor required to account for the positron and/or electron plus positron excess via annihilation to 
muons is compatible with present BBN constraints. 

We next consider constraints from gamma-rays. A possible gamma-ray signal from the GC due 
to DM annihilation has been extensively studied as it could provide a good method to study the 
nature of DM astrophysically. So far the HESS group has reported that power-law signals were 



observed from the GC 1145 . 



46j] for 200 GeV < E y < 700 GeV. Quite recently the FERMI satellite 



group also reported their preliminary result for the signals observed from the galactic mid-latitude 
(10° < \b\ < 20°) for 200 MeV < E y < 10 GeV. When we adopt a cuspy profile of the galaxy, such 
as the NFW profile, the gamma-ray signal from muon emission can exceed the observed signal. 
However, if we take a milder profile such as the cored isothermal profile, then for the moment DM 
annihilation is not constrained by the current observations [3(1 31, 47] 9 . To clarify the dependence 
of the DM constraints on the density profile, we need more accurate data on the diffuse gamma-ray 
background, which will be provided by FERMI in the near future. 

In addition, there are CMB constraints on the enhanced S'S annihilation cross-section. It has 
been shown in ref. I49I that energetic particles from rapid S^S annihilation can reionize neu- 



tral hydrogen at the last scattering surface, leaving an imprint on the CMB. The analysis of [48] 
concludes that current data from WMAP5 imposes a 2-o upper bound on the S^S annihilation 
cross-section which is given by 

3.6 x 10- 24 cm 3 /s / M s \ 



<° v W x t X3 < y ^— ) , (36) 



9 Note that the positron and electron plus positron signals will not change even if we used the cored isothermal 
profile because local annihilation within 1 kpc dominates the production of electrons and positrons with > 10 GeV 
energies. 
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where / is in the range 0.2 — 0.3 for annihilation to 2 /u + /u pairs. Thus a boost factor of 0(1000) 
is marginally allowed by the current data. 

Finally, we briefly comment of the possibility of neutrino signals from the GC. Detecting such 
neutrino signals in the future might be useful to distinguish the Higgs portal DM model from 
others, since muon neutrinos are produced by the decay of the pairs coming from DM 

annihilation and subsequent %' 3 decay. So far Super-K has reported upper bounds on the up-going 
muon flux coming from neutrinos emitted from the GC Q50U . We can compare the theoretical 
prediction of the neutrino flux in our model with this Super-K upper bound. According to the 



discussion of Ref. plQ . our model is presently allowed since neutrinos are not produced directly 



but indirectly through the decay of the charged le pton s and possibly mesons. It is expected that 



future neutrino experiments such as KJVBNeT 11521 



53Q or IceCube DeepCore 11541 



55Q will be able 



to detect the up-going muons induced by the neutrinos emitted from the GC. 



III. A SOMMERFELD ENHANCED VERSION OF THE MODEL 

In the previous section we studied the conditions for a successful Higgs portal model with 
non-thermal production of DM. In this section we will replace non-thermal production with Som- 
merfeld enhancement of the annihilation cross-section as the source of the boost factor. The main 
difference between the two models is the reduced number of hidden sector scalar fields, since 
%iis no longer needed as the source of the non-thermal DM density. This will also eliminate the 
most heavily suppressed O(10~ 10 ) couplings, which were necessary to ensure the late decay of %\. 
The S DM annihilation to %3 pairs and subsequent %3 decay to pairs is unchanged from the 
non-thermal scenario. 

Since in the non-thermal model there must exist a light scalar %3 if we wish to avoid leptophilic 
couplings, it is natural to ask whether we can eliminate %i and consider instead thermal DM 
with a Sommerfeld enhanced annihilation cross-section, with the attractive force mediated by %3- 
exchange. The correct thermal relic density of S DM is obtained if the S^Sy^Xi coupling is in 
the range 0.1-1 for Ms ~ 0.1 — 1 TeV J24 ]. If we then consider the coupling (^5)23 in Eq.(l) and 
introduce < >, we obtain the interaction 

(r|s) 23 <%2> XsS^S + h.c. . (37) 

This interaction can produce the required long-range force between S particles via %3 exchange. 
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The condition for a Sommerfeld enhanced annihilation rate is M X3 < CcMs, where a = X 2 /4n and 
the effective coupling from %3 exchange is A, « (^5)23 < %2 > /M$. Therefore 

Since M %3 ~ 200 MeV in our model, this will be satisfied if (^5)23 > 0.4 when M s ~ 1 TeV. 

Therefore, in addition to simplifying the model by eliminating %i, Sommerfeld enhancement 
permits larger DM masses, Ms ~ 1 TeV. This may be significant in light of recent analyses [28, 29I 
of the new FERMI and HESS electron plus positron data, which favour DM particles with TeV 
scale masses which annihilate to muons (with the case of annihilation to 4/j being favoured by 
the analysis of D29J]). Since the Higgs portal models generally predict that DM annihilates to two 
pairs via decay of the primary %3 pair, a Sommerfeld enhanced version of the Higgs portal 
model, in contrast with the non-thermal model, could provide an explanation for both the higher 
energy electron plus positron excess and the lower energy PAMELA positron excess. 



IV. CONCLUSIONS 



We have considered two DM models for cosmic ray excesses which are based on Higgs portal- 
type couplings of a scalar DM sector to the SM, one with non-thermal DM as the explanation of 
the boost factor and the other with thermal DM and Sommerfeld enhancement of the annihilation 
cross-section. 

In the case of the model with non-thermal production of DM, the DM scalar mass must be less 
than about 600 GeV if the model is to remain perturbative. Therefore if this model is correct then 
the PAMELA positron excess can be explained by DM annihilation but the higher energy electron 
plus positron flux suggested by FERMI and HESS must have a different explanation. This is likely 
to be true of most models without Sommerfeld enhancement. Non-thermal production of DM is 
possible via quartic scalar couplings. However, the couplings leading to decay of the heavy scalar 
which produces the DM density must be highly suppressed, < 10~ 10 , in order to ensure that the 
heavy particle decays well after the DM particle freezes-out. 

A successful model must also account for DM annihilation to primarily leptonic states. If we 
do not wish to introduce DM which couples preferentially to leptons then the only way to achieve 
this is kinematically, by ensuring that DM annihilates to unstable final states which are too light 
to subsequently decay to hadrons. Our model can generate the required decay process via mixing 



19 



of the %3 scalar of the hidden sector with the Higgs, leading to the decay of %3 primarily to n + iT 
via the muon Yukawa coupling if its mass is in the range 2m M — 2m % o (212 — 270 MeV). The small 
%3 mass requires that the quartic scalar couplings of %3 to the Higgs and to %2 are ~ 10~ 6 . The 
final state is essential if we require that the %3 density decays prior to nucleosynthesis (which 
%3 would otherwise dominate) but does not decay to pions or nucleons, which would produce a 
potentially dangerous photon or antiproton flux. This is a clear prediction of the Higgs portal 
model, which applies equally to the Sommerfeld enhanced version. 

We conclude that quartic couplings of a relatively simple scalar DM sector can achieve the 
required enhancement of the annihilation rate and leptonic final states, but appropriate mixtures 
of strongly suppressed and unsuppressed quartic couplings and large and small mass terms are 
required. In the absence of symmetries or dynamical effects which can explain them, such hierar- 
chies would appear unnatural. It is therefore to be hoped that the pattern of masses and couplings 
can be understood in terms of the symmetries or dynamics of a complete theory, for which the 
present model is the low energy effective theory. 

A significant feature of the Higgs portal model, which can distinguish it from those with purely 
leptophilic annihilation modes, is that there can be a significant coupling of DM to Higgs pairs. 
This could produce a non-negligible antiproton component in the cosmic rays from DM annihila- 
tion if the annihilation process S^S — > H^H is not too suppressed relative to the dominant process 
->■ %3%3. The S^SH^H coupling may also allow direct detection of DM. 

Constraints from BBN are important for the model with non-thermal DM, since the annihilation 
rate is large at all temperatures. We found that both the muon and Higgs final states are consistent 
with an annihilation cross section as large as 10~ 23 cm 3 s _1 for Ms ^ 600 GeV. The model is also 
consistent with the gamma-ray signal from the galactic centre and from the diffuse gamma-ray 
background in the case of a cored isothermal halo profile, but not in the case of a cuspy NFW 
profile. 

In the Sommerfeld enhanced version of the model, the low mass scalar which accounts for 
leptonic DM annihilation also mediates the force responsible for the Sommerfeld enhancement. 
In this case we can reduce the number of additional scalars by one, since %i is no longer needed to 
produce the DM non-thermally, which also eliminates the most highly suppressed couplings. This 
version of the model can accomodate a TeV scale DM particle, allowing it to explain the electron 
plus positron excess suggested by FERMI and HESS as well as the positron excess observed by 
PAMELA. Exactly as in the non-thermal model, DM annihilates to %3 pairs which subsequently 
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decay via Higgs mixing to 2 fj + fj pairs. This may be significant, as recent analyses suggest that 
the new FERMI and HESS electron plus positron data favours TeV scale DM particles annihilating 



to muons [|28 



L29] 



290, with annihilation to 4/u via intermediate decaying scalars being favoured by 



The Higgs portal models considered here should have phenomenological signals due to the 
coupling of the DM sector to the Higgs bilinear and the mixing of the Higgs with the SM singlet 



%2 Q56(] . If S or %2 are light enough then they may be produced via Higgs decay at the LHC. 
The mass eigenstate Higgs boson is also be expected to have a large singlet component, with 
consequences for Higgs phenomenology. These features may not be unique to our model, but they 
would provide indirect support for it. In addition, the muon neutrinos produced by the decay of the 
pair from DM annihilation may be detectable via upward-moving muons at future neutrino 
experiments. 

Acknowledgement 

The research at Lancaster is supported by EU grant MRTN-CT- 2004-503369 and the Marie 
Curie Research and Training Network "UniverseNet" (MRTN-CT-2006-035863). 



[ 1 ] O. Adriani et al. , larXiv:0810.49"95l [astro-ph] . 

[2] S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482, L191 (1997) 
| arXiv:astro-ph/9703192) ; J. J. Beatty et al, Phys. Rev. Lett. 93 (2004) 241102 
| arXiv:astro-ph/0412230) . 

[3] M. Aguilar et al. [AMS-01 Collaboration], Phys. Lett. B 646, 145 (2007) | |arXiv:astro-ph/0703 154) . 

[4] O. Adriani et al. , larXiv:08 10.49941 [astro-ph] . 

[5] J. Chang et al, Nature 456, 362 (2008) . 

[6] S . Torii et al. , larXiv:0809.0760l [astro-ph] . 

[7] F. Collaboration, larXiv:0905.0025l [astro-ph.HE] . 

[8] H. E. S. Aharonian. [arXlv:0905.0105i rastro-ph.HE1. 

[9] D. Hooper, A. Stebbins and K. M. Zurek. larXiv:0812.3202 [hep-ph]. 



21 



[10] D. Hooper, P. Blasi and P. D. Serpico, JCAP 0901, 025 (2009) la rXiv:0810. 15271 Tastro-phll ; S. Pro- 
fumo, arXiv:08 12.44571 rastro-phl. 

[11] H. Yuksel, M. D. Kistler, and T. Stanev. larXiv:08 10.27841 K. Ioka. larXiv:0812.4851 ; N. J. Shaviv, 
E. Nakar, and T. Piran, larXiv:0902.0376l E. Borriello, A. Cuoco and G. Miele, arXiv:0903. 18521 
[astro-ph.GA]; P. Blasi. larXiv:0903.2794l N. Kawanaka, K. Ioka and M. M. Noiiri. larXiv:0903.3782l 
[astro-ph.HE]; Y. Fujita, K. Kohri, R. Yamazaki and K. Ioka. larXiv:0903.5298 1 rastro-ph.HE] ; P. Blasi 
and P. D. Serpico, arXiv:0904.087T1 

[12] S. Shirai, F. Takahashi and T. T. Yanagida. larXiv:0905.0 388 [hep-ph]; C. H. Chen, C. Q. Geng and 
D. V. Zhuridov. larXiv:0905.0652l rhep-phl . 

[13] Y Nomura and J. Thaler. larXiv:08 10.53971 rhep-phl . 

[14] M. Cirelli, M. Kadastik, M. Raidal and A. Strumia. larXiv:0809.2409l rhep-phl . 
[15] J. Hisano, S. Matsumoto, O. Saito and M. Senami, Phys. Rev. D 73, 055004 (2006) 
| |arXiv:hep-ph/0511118) . 

[16] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer and N. Weiner, Phys. Rev. D 79, 015014 (2009) 



[arXiv:08 10.0713 [hep-ph]]. 
[17] M. Pospelov and A. Ritz, Phys. Lett. B 671 (2009) 391 HarXiv:0810.1502i rhep-ph11. 



[18] T. Moroi and L. Randall, Nucl. Phys. B 570, 455 (2000) [arXiv:hep-ph/9906527 ]; M. Fairbairn and 

J. Zupan. larXlv:0810.4147 rhep-phl . 
[19] J. McDonald. larXiv:0904.0969l rhep-phl . 

[20] K. Enqvist and J. McDonald, Phys. Lett. B 440, 59 (1998) flarXiv:hep-ph/9807269) ; M. Fujii 
and K. Hamaguchi, Phys. Lett. B 525, 143 (2002) | |arXiv:hep-ph/01 10072] ; Phys. Rev. D 66, 
083501 (2002) larXiv:hep-ph/0205044) ; M. Fujii and M. Ibe, Phys. Rev. D 69, 035006 (2004) 
[arXiv:hep-ph/0308118]. 



[21] D. Feldman, Z. Liu and P. Nath, Phys. Rev. D 79, 063509 (2009) [arXiv: 08 105762 [hep-ph]]. 

[22] M. Ibe, H. Murayama and T. T. Yanagida, Phys. Rev. D 79, 095009 (2009) I arXiv:08 12.0072 Qiep- 



ph]]; W. L. Guo and Y. L. Wu, Phys. Rev. D 79, 055012 (2009) [arXiv:090 1.1450 [hep-ph]]; M. Ibe, 
Y. Nakayama, H. Murayama and T. T. Yanagida, JHEP 0904, 087 (2009) HarXiv:0902.2914i rhep-ph11. 



F. Y. Cyr-Racine, S. Profumo and K. Sigurdson, arXiv:0904.3933 [astro-ph.CO]. 
[23] P. J. Fox and E. Poppitz, larXiv:08 11.03991 [hep-ph]. H. S. Goh, L. J. Hall and P. Kumar, 

larXiv:0902.08l4l [hep-ph] . 
[24] J. McDonald, Phys. Rev. D 50, 3637 (1994). 



22 



[25] V. Silveira and A. Zee, Phys. Lett. B 161, 136 (1985); C. P. Burgess, M. Pospelov and T. ter Veldhuis, 
Nucl. Phys. B 619, 709 (2001) | |arXiv:hep-ph/0011335) . 



[26] B. Patt and F. Wilczek, arXiv:hep-ph/0605188 

[27] M. C. Bento, O. Bertolami, R. Rosenfeld and L. Teodoro, Phys. Rev. D 62, 041302 (2000) 



[arXiv:astro-ph/0003350|; M. C. Bento, O. Bertolami and R. Rosenfeld, Phys. Lett. B 518, 276 



(2001) [arXiv:hep-ph/0103340 ]; J. McDonald, Phys. Rev. Lett. 88, 091304 (2002); L. Lopez Honorez, 
E. Nezri, J. F. Oliver and M. H. G. Tytgat, JCAP 0702, 028 (2007) | |arXiv:hep-ph/06 1 2275 ) ; N. Sahu 
and U. Sarkar, Phys. Rev. D 76, 045014 (2007) larXiv:hep-ph/0701062 ] ; J. McDonald, N. Sahu 
and U. Sarkar, JCAP 0804, 037 (2008) HarXiv:071 1.48201 [hep-phll; J. March-Russell, S. M. West, 
D. Cumberbatch and D. Hooper, JHEP 0807, 058 (2008) HarXiv:0801.3440l rhep-phll; J. McDon- 
ald and N. Sahu, JCAP 0806, 026 (2008) l arXiv:0802.3"847l [hep-ph]] ; Y. G. Kim, K. Y. Lee and 
S. Shin, JHEP 0805 (2008) 100 HarXiv:0803.2932l rhep-phll ; D. G. Cerdeno, C. Munoz and O. Seto, 
larXiv:0807.3029l [hep-ph] . 
[28] L. Bergstrom, J. Edsjo and G. Zaharijas, arXiv:0905.0333 [astro-ph.HE]. 



[29] P. Meade, M. Papucci, A. Strumia and T. Volansky, arXiv:0905.0480 [hep-ph]. 



[30] G. Bertone, M. Cirelli, A. Strumia and M. Taoso, arXiv:08 11.3744 [astro-ph]. 
[31] M. Kawasaki, K. Kohri and K. Nakavama. larXiv:0904.3626l [astro-ph.CO1. 
[32] K. Griest and M. Kamionkowski, Phys. Rev. Lett. 64, 615 (1990). 
[33] B. W. Lee and S. Weinberg, Phys. Rev. Lett. 39, 165 (1977). 

[34] T. Delahaye, R. Lineros, F. Donato, N. Fornengo and P. Salati, Phys. Rev. D 77, 063527 (2008) 
HarXiv:0712.2312l [astro-ph]] . 

[35] M. Cirelli, R. Franceschini and A. Strumia, Nucl. Phys. B 800, 204 (2008) HarXiv:0802.33781 [hep- 
ph]]. 

[36] I. V. Moskalenko and A. W. Strong, Astrophys. J. 493, 694 (1998) | |arXiv:astro-ph/9710124) . 
[37] E. A. Baltz and J. Edsjo, Phys. Rev. D 59, 02351 1 (1999) | |arXiv: astro-ph/9808243 ) . 
[38] J.F Navarro, C.S. Frenk and S.D.M. White, Astrophys. J.462, 563 (1996) . 

[39] P. Gondolo, J. Edsjo, L. Bergstrom, P. Ullio, M. Schelke and E.A. Baltz, |astro-ph/0406204 ]. Dark- 



SUSY homepage http://www.physto.se/~edsjo/darksusy/JCAP 0407 (2004) 008. 
[40] K. Jedamzik, Phys. Rev. D 70, 083510 (2004). 

[41] J. Hisano, M. Kawasaki, K. Kohri and K. Nakavama. larXiv:08 1 0. 1 8921 [hep-phi . 

[42] J. Hisano, M. Kawasaki, K. Kohri, T. Moroi and K. Nakavama, rarXiv:0901 .358 2 [hep-ph]. 



23 



[43] M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 625, 7 (2005); Phys. Rev. D 71, 083502 (2005). 
[44] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006). 
[45] F. Aharonian et al. [The HESS Collaboration], Astron. Astrophys. 425, L13 (2004); 
[46] F. Aharonian et al. [H.E.S.S. Collaboration], Phys. Rev. Lett. 97, 221102 (2006) [Erratum-ibid. 97, 
249901 (2006)]. 

[47] P. Salucci, A. Lapi, C. Tonini, G. Gentile, I. Yegorova and U. Klein, Mon. Not. Roy. Astron. Soc. 378, 
41 (2007) |arXiv:astro-ph/0703 115) ; M. Persic, P. Salucci and F Stel, Mon. Not. Roy. Astron. Soc. 
281, 27 (1996) [ |arXiv:astro-ph/9506004) . 

[48] T. R. Slatyer, N. Padmanabhan and D. P. Finkbeiner, Phys. Rev. D 80, 043526 (2009) 
HarXiv:0906. 11971 [astro-ph.CO]] . 

[49] T. Kanzaki, M. Kawasaki and K. Nakavama. irXiv:0907.3985l rastro-ph.COl. 

[50] S. Desai et al. [Super-Kamiokande Collaboration], Phys. Rev. D 70, 083523 (2004) [Erratum-ibid. D 

70, 109901 (2004)]; Astropart. Phys. 29, 42 (2008). 
[51] J. Hisano, M. Kawasaki, K. Kohri and K. Nakavama. larXiv:08 1 2.02191 [hep-phi . 
[52] Kappes, A. et al. 2007, ApJ, 656, 870 
[53] Kappes, A. & Consortium, f. t. K. 2007. larXiv:071 1.0563 1 
[54] Achterberg, A. et al. 2007, PRD, 76, 042008 [Erratum: PRD, 77, 089904 (2008)] 
[55] D. F. Cowen [IceCube Collaboration], J. Phys. Conf. Ser. 110, 062005 (2008). 

[56] V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf and G. Shaughnessy, Phys. Rev. D 77, 
035005 (2008) liarXiv: 0706 .431 1 1 rhep-phll. 



24 



